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Abstract 

We extend our studies of a quantum field model defined on a lattice having the dilation 
group as a local gauge symmetry. The model is relevant in the cross-disciplinary area 
of econophysics. A corresponding proposal by Ilinski aimed at gauge modeling in non- 
equilibrium pricing is realized as a numerical simulation of the one-asset version. The 
gauge field background enforces minimal arbitrage, yet allows for statistical fiuctuations. 
The new feature added to the model is an updating prescription for the simulation that 
drives the model market into a self-organized critical state. Taking advantage of some 
flexibility of the updating prescription, stylized features and dynamical behaviors of real- 
world markets are reproduced in some detail. 
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1. Introduction 

The analysis and modeling of financial price time series has a long history [H, [3j S| 
and has attracted considerable interest at an accelerated pace in the last two decades. 
Technological advances have made it possible to collect and process vast amounts of 
data. As a result, various stylized facts about the statistics of financial data have been 
discovered 0, [^, @|- These features are mostly concerned with scaling laws, akin to 
findings in many systems described by statistical physics. There, scaling behavior arises 
from the interaction of many units in such a way that a critical state is reached. Thus, one 
may ask if a financial market, for example, can be modeled based on similar principles. In 
generic terms, the building blocks could be many individual agents with suitable mutual 
interactions. Indeed, Lux and Marchesi have shown that scaling laws can arise in such a 
setting 0. 

When building a microscopic model it is prudent to rely on a theoretical foundation 
supported by evidence. In the present work, we will employ two such principles. First, ar- 
bitrage opportunities will be annihilated during the time evolution of the market, though 
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admitting statistical fluctuations. Second, the dynamics of the model will drive it into a 
self-organized critical state, thus naturally giving rise to scaling behavior. Both aspects 
have been investigated separately in previous work, see |8] and 9] respectively. Here, we 
merge those elements into a microscopic market model, using numerical simulation to 
study its characteristics. 

The next section gives an overview of the model's dynamics and definitions. 

2. Lattice model 

Following a proposal by Ilinski [l^] we define a lattice field theory with a ladder 
topology as shown in Fig.[TJ In physics terms there are matter fields ^{x) £ R+ defined 
on sites x = where j — . . .n means discretized time, and i = 0, 1 is a spatial 

index. These are represented by filled circles in Fig. [TJ In addition, there are gauge fields 
Ofj,{x) G K"*" which live on links starting at site x in temporal fi — or spatial /i = 1 
direction. Those are represented by arrowed lines in Fig. [T] 




Figure 1: Left: Illustration of the ladder geometry of the lattice model and the label scheme for the 
fields. Right: Depiction of the gauge invariant elements A, B, C, D used in the action. 

As a model for a financial market, again following [lo| . we interpret the matter field as 
instances of an account value, in some unit. At « = it could be a cash account, whereas 
at i = 1 the value may be interpreted as the number of shares owned in some financial 
instrument. The spatial links Qi(0,j), connecting cash and shares, are simply conversion 
factors between the corresponding units. Temporal links 0o(O,j'), which connect cash 
value sites one time step apart, mean interest rate factors. Similarly, temporal links 
0o(l,j), starting from a shares site, carry information about the change in share value 
one time step apart. 

The rationale behind such a model is to describe a market that dynamically evolves 
independent of the trading units being used. For example, in comparable markets, the 
dynamics should not depend on the specific, notably arbitrary, currency unit being used 
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in transactions. This, at least, is the hypothesis which should apply to markets trading 
in like instruments. 

Mathematically, this idea is implemented by a quantum field theory with local gauge 
invariance. Such has been worked out in great detail in the context of financial markets 
[lo| . In a previous work we have studied some aspects of those ideas using numerical 
simulation J8| . Since the current work is directly building on the latter, we refer the 
reader to [8| for the technical details. In particular, we shall use the nomenclature 
therein. However, to keep the presentation self contained, the essential building blocks 
are discussed in what follows. 

The dynamics of the model derives from an action 5[6, $, i>] for the lattice fields that 
is invariant with respect to local gauge transformations 

^ gixMx) (1) 
^x) ^ Hx)g-\x) (2) 
e^(x) ^ g{x)&^{x)g-^ix + ef,), (3) 

where $(a;) — l/$(a;) and g{x) £ G is an element of the dilation group G — K+, i.e. 
multiplication by positive real numbers. Those carry out conversions between (arbitrary) 
units. The action is constructed from the elements depicted in Fig. [TJ These are the 
smallest gauge invariant objects that can be assembled from the fields. 

The diagram associated with label A is known as the elementary plaquette 

P^^Ax) - e^(x)e,(a; + e^)Q-'{x + e,)e-\x) . (4) 

Its value is interpreted as the gain (or loss) realized through an arbitrage move |10|. 
The global minimum of the classical action 5[8, $, $] corresponds to zero arbitrage [8|. 
Quantization of the field is done through the usual path integral formalism. The partition 
function thus is defined as the functional integral 

Z(/3) =y"[i:>e][L>$]e"/5'^[®'^'^] . (5) 

In this way stochastic fluctuations about zero arbitrage are allowed. Their magnitude is 
regulated by the parameter /3. 

Diagrams B,C,D are gauge invariant elements of the form 

Rf,{x) =^x)e^{x)<^{x + ef,), (6) 

where is a unit vector in direction /i. Diagram C, for example, gives the value of the 
investment instrument at time j + 1 divided by its value at j, provided we adopt the 
above interpretation of the fields. It is a measure for the relative change of the asset 
value during one time step 

Ro{l,j) = l>(l,j)eo(l,j)$(l,j + 1) . (7) 

We also define the related quantity 

rj+i -logi?o(l,j), (8) 

commonly called the 'return', indicating a gain (> 0) or a loss (< 0) at the end of the 
time step. 
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3. Updating strategy 



The generation of lattice field configurations as implemented in [8| follows a standard 
procedure. Based on the action S'[9,<i>,l>] the field components are updated through 



a heatbath algorithm [111 . Il2j linked to the partition function ([S]). Periodic boundary 
conditions (in the time direction) are imposed on all fields as well. However, in contrast 
to 0, the updating strategy is modified in two respects. 

First, we do set constraints on the fields that live on the axis i — 0, see Fig. [TJ The 
reasoning here is that we wish to design the model such that the axis describes a cash 
account subject to accumulating interest. The interest rate is endogenously determined. 
Even at 10% annually the daily rate factor is 1.0003 and thus hardly distinguishable 
from one. As the model is designed to describe a high frequency market, where the time 
extent n translates to typically a day, or so, we wish to set a constraint accordingly. In 
a gauge model this is not straightforward, because the meaning of the field components 
is gauge dependent. To remedy this situation, gauge fixing is called for. With reference 
to (IT][3]) define a gauge transformation along the axis i — 0, 

5(0,j) = $(0,j), (9) 

with g{x) on all other sites being arbitrary. The gauge transformed fields along the axis, 
i = 0, then are 

^'{0,J) = 5(0,j)<i>(0,i) = l (10) 
$'(0,j) = $(0,j)g-i(0,j) = l (11) 

e'o{o,j) - $(o,j)eo(o,j)$(o,j + i) = i?o(o,j). (12) 

In the last equation we recognize the link variable as the (gauge invariant) return of 
the cash holding during one time step. We therefore set 

e'o{o,j) = i. (13) 

In our simulation we choose a random start for the lattice fields. From there, the con- 
straint i?o(0,j) = 1 is then implemented by applying the gauge transformation ([9|), and 
then setting the axis links to one (|13p . During the subsequent updating procedure the 
axis fields $(0, j) and 8o(0,j) are never changed. Nonetheless, the right-hand side of 
(I13p may differ from one, depending on the interest rate factor desired. 

The next step is to run a heatbath algorithm with the lattice action ^[B, $, $] until 
equilibrium is reached Q. The lattice field configurations then model a market environ- 
ment where arbitrage opportunities exist only briefly, subject to fluctuations due to the 
quantum nature of the fields. In economic terms this model describes an efficient market. 
Equilibrium, however, does not seem to be realized in the real world [13] . 

Second, subscribing to this paradigm, we introduce a new element which is applied 
post equilibrium. In (9| we have studied a simple model where market instances live 
along a linear chain in time direction. The sites carry fields rj S which are directly 
interpreted as returns, thus having the same meaning as There is no gauge field in 
the simple model. The key ingredient is an updating algorithm that mimics the popular 
Bak Sneppen evolutionary model [l4, 15, 1^. The quantity 
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(14) 



turned out to be essential to the field dynamics. In the context of [9*] the updating 
strategy consists in finding the site js for which the absolute value of ([T4| is maximal 
\vj^ I = maxjlvjl : j = . . . n}, and then replace rj^ and its neighbors rj^±i with random 
numbers. Iterating this process leads to a self-organized critical state and produces price 
times series, and related statistics, which are almost indistinguishable from those in real 
markets 01 . 

In view of those results it seems desirable to replicate this updating strategy within 
the framework of the gauge model as closely as possible. Towards this end, we still do 
define the 'fitness' measure Vj as in ([HI) , however, the returns are now given by ([5]), and 
([7]) . Their composition is illustrated in Fig. [51 The updating prescription proceeds with 
finding the 'signal' 

V — max{|wj| : j = 1 . . .n} (15) 
of the field configuration, and the site js of its location 



= V. 



(16) 



l),eo(l,j. - l),$(l,Js), which 
J, , and the two next-neighbor links <do{l,js),Qo{l,js — 2), see Fig. [21 
+ 1) and ^{l,js — 2) are left unchanged. In this way only the three 



We then update the three field components $(l,j, 
enter the return r 
Note that $(l,j. 



returns rj_, _ i , rj^ , rj^ 
prescription used in 



are affected. This strategy most closely resembles the updating 
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Figure 2: Left: Illustration of the returns involved in the 'fitness' criterion II14I I (left), and the field 
components subject to updating (right), done at the 'signal' site j = js- 



Updating those field components consists in drawing random numbers from certain 
probability distributions. We have chosen those based on the lattice action 5'[6, $, $] 
mentioned above. Heatbath steps using the corresponding Boltzmann-like distribution, 
see ([S]), are applied to the various field components. Essentially, the probability distri- 
bution for a given field component is given by its local environment. It is convenient to 
rewrite the fields as 

e,,(x) = e^f'(^) and $(x) = e'^(^) . (17) 

Then, after a gauge transformation, the probability densities for the gauge fields and the 
matter fields, respectively, have the form 

Peie^ix)) cx exp(-2/3x/i0Le cosh(0^(x)) ) (18) 
p$(0(a;)) oc exp(— 2/3-\/ i$Z$ cosh(0(a;)) ) . (19) 
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These results are derived in detail in Appendix [Appendix A[ The coefficients Lq,Lq 
and L^,Lq, are independent of 6^(a;) and <&(a;), respectively. The products LqLq and 
L$L$ are gauge invariant and, together with the parameter ^, determine the variance 
of the probability distributions for the field components. Those distributions strongly 
depend on the local environment at the location of the fields. 

Now, at each updating step we randomly draw fields (^'(1, j), js — 1 < j < js, from 
(IT9)) . Relevant averages considered are 

js js 

ae = l J2 and = ^ £ (20) 

j=js-'2 j=js-l 

Updating the fields then is accomplished by replacing 

0o(l,j) ^ 0o{l,j)~xad, j.'2<j<j, (21) 
0(1, j) ^ 0'(l,j)-a0, Js ~l< J <js. (22) 

By including the parameter x we have introduced a novel feature to the updating process. 
While X = 1 essentially mirrors the strategy in 9] , deviations from that value introduce 
very interesting features to the model. We will be able to describe a range of different 
returns distributions and time series, as will be described in the next section. 

Finally, we apply heatbath updates to the two spatial (horizontal) link variables 
0i(O, j), js — 1 < j < js, which connect to the affected matter fields, see (P^. These 
links occur in three elementary plaquettes 

Pio(o,i - 1) = ei(o,j - i)eo(i,i - i)e^\o,j)eo\o,j - 1) , (23) 

see ([1]), where — 1 < j < + 1- The reason is that updating 0o(l, j), as prescribed by 
(PT|) . changes the plaquettes ([25)) and thus upsets the no-arbitrage environment of the 
lattice fields. Updating the above links with the lattice action S[Q, <i>, $] rectifies this 
circumstance. 

4. Results 

The simulations were done on a lattice of size n = 782 with gauge field coupling 
parameter (3 = 1, and the matter field couplings = = 1. These parameters are the 
same as in Q, with the one asset model m — 1. The number of heatbath update steps 
was 10^ to equilibrate the field from a random start. Final configurations were reached 
after 4 x 10^ 'signal' updates. 

First, we discuss the effect of the parameter x in pl]l . A suitable observable (order 
parameter) is the gauge invariant link along the asset axis ([7]). Using the notation pT|) 
we have 

i?o(l, J - 1) = exp(-(/.(l, J - 1) + 0o(l, J - 1) + 0(1, j)) = exp(r,) . (24) 

The updating algorithm, described in Sect. [31 employs symmetric probability distribution 
functions for 9 and 4>. Consequently, the probabilities for realizing a gain rj > and a 
loss rj < are equal. Thus we define the symmetric link 

= I (exp(rj) + exp(-rj)) - 1 = cosh(rj) - 1 (25) 
^ 6 



and its lattice average 

1 " 

L=-Y.L,. (26) 
i=i 

Numerically, the value of L is particular to a distinct lattice field configuration. We denote 
the (stochastic) average over field configurations with angle brackets, here {L). In Fig. [3] 
the dependence of {L) on the parameter x is displayed. The plot symbols '•' indicate 
data points from simulations at x = 10*^, k = —6, —5 . . . + 1. The errors come from 48 
field configurations. The line curve in Fig. |3]is a four parameter fit, ai . . .04, to those 
eight data points with y = ai tanh[a2(— log]^g(a;) + 03)] + 04. Remarkably, there clearly 
is a transition region in the, approximate, range 10~* < X < 10^^. For small values 
X <C 10^'' the average link operator saturates at 0.2406, while for large values 10^^ <C x 
it tends to 0.0092, according to the fit. Using the (crude) conversion cosh(r) — 1 = L for 
simplicity, see (pS)) . this corresponds to returns r sa ±0.68 and r « ±0.14, respectively, 
for the above limits of x- Those limits may describe valid markets, which could be seen as 
volatile and calm, respectively. In view of this observation the transition region becomes 
particularly interesting. It opens up the possibility to simulate markets with a wide range 
of features between those extremes. Below, we will present results for x — 0.0005, 0.0013 
within the transition region. Those are indicated by plot symbols '+' in Fig. [3) 
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Figure 3: Expectation value (L) of the gauge invariant average link operator II26I I as a function of the 
update parameter x in I I21II . The plotting symbols '•' and '+' correspond to specific values of x for 
which simulations were made. The symbol 'o' indicates the symmetry point x = 0.0052 of the fit. 

During a simulation, the evolution of the lattice towards a critical state can be mon- 
itored, for example, by observing the signal V, see (ITSt . as a function of the updating 
step counter, say s = 0, 1, 2 . . .. Writing V{s) we follow 16] and define the 'gap' function 

G{x) = min{V^(s) : s e N U {0} and s < x} with a; £ K+ U {0} . (27) 

This is a decreasing piecewise constant function with discontinuities at certain discrete 
values Xfc, fc G N. The set of update steps between Xk-i and Xk is called an avalanche 
of length Ak = Xk — Xk-i- Eventually, as x — >■ cx), the avalanche size diverges and the 
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system has reached criticahty [16|. For an elaboration on these concepts, presented in a 
context close to the current work, we refer the reader to We here only show a key 
result. 

In Fig. U the frequency distribution AiV/AA of the avalanche sizes is displayed. Here 
AA is a binning interval for the avalanche sizes and AA^ is the count of avalanches within 
that interval. We have used 10000 bins with A A = 1. The data points and errors come 
from an ensemble average over 2400 independent lattice simulations with 4 x 10® update 
steps each. The log-log plot clearly shows power law behavior. A power law indicates 
scaling, which is a signature feature of a critical system. 




Figure 4: Frequency distributions of avalanche sizes for two values x = 10 *,10 of the update 
parameter. 

For the time being, we continue to present results for the two update parameters 
X = 10~^ and % = 10""'^. These values correspond to the boundaries of the transition 
region, see Fig. [31 Inside that region, the frequency distributions of avalanche sizes are 
almost indistinguishable from the results shown in Fig. |4l Examples of model markets 
for suitable update parameters in the transition region are discussed below. 

The gains distributions produced by the lattice model with x — lO^'^, 10^^ are shown 
in Fig. [S] The gauge invariant returns r, as defined in ([5]), are put into bins of size 
Ar, and Ac/Ar is the number of counts per bin. The errors are obtained from 2400 
independent simulations. While both histograms possess fat tails, we observe a distinct 
difference of the qualitative features for the distributions. At x = 10^^ a distinctly 
pointed central peak sits on very broad bulging tails. Looking at x — 10^^ the central 
peak has broadened such that its very top is almost Gaussian while the tails look linear. 
This means that the two distributions cannot be mapped into each other by simple 
scale transformations applied to the axes. The two gains distributions describe genuinely 
different markets. Interestingly, this matches our previous assessment of the the regions 
X ^ 10~^ and 10~^ ^ x discussed in the context of Fig. [31 

Examples of returns time series for each of the two parameters x ^-re displayed in 
Fig. [HI By visual inspection, both of those exhibit volatility clustering, yet display 
different dynamical behavior. In the realm of % = 10~^ the time series appears to favor 
one side of the zero mark for short periods of time, as compared to the series of x = 10~^ 
which smoothly fluctuates about zero in either direction. 

For a closer investigation, we have selected the parameters x = 0.0005, 0.0013 from 
the transition region. In Fig. [31 their locations are marked by '+' plot symbols. The 
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Figure 5: Lattice generated gains distributions of the returns r for the values x = 10 ^,10 ^ of the 
update parameter. 
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J J 
Figure 6: Examples of returns time series rj versus the lattice time j at x = and x = 10~^. 

choice of these values reflects our observation that model market characteristics, say the 
gains distribution for example, hardly change as x is decreased from « 0.1 to « 0.005. 
Most of the market model tuning happens in the upper segment of the (L) curve in Fig. [3] 
as X further decreases below « 0.005. 

It is tempting to utilize this range for modeling historical markets of varying char- 
acteristics. However, we leave this to future work. Here, our focus is on the properties 
of time series dynamics of the lattice model as being evaluated with standard financial 
analysis tools. 

The gains distributions for x = 0.0005 and x = 0.0013 are displayed in Fig. [T] 
They were obtained from 2400 simulation runs with length 4 x 10^ each. Again, the 
distributions clearly exhibit fat tails but are otherwise different in their shapes. The 
X = 0.0005 distribution puts more emphasis on larger returns than the x = 0.0005 
distribution. This is echoed in the corresponding returns time series. 

Samples of those time series are displayed in Fig. |S] for x = 0.0005 and Fig. [3] for 
X = 0.0013, respectively. Each figure is composed of eight randomly selected independent 
simulation runs with 4 x 10^ updates each. Compared to the x = 0.0013 time series, the 
X — 0.0005 series have a, somewhat perceptible, higher occurrence of volatility clusters 
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Figure 7: Gains distributions from simulations with update parameter x = 0.0005 and x = 0.0013. 

as well as bigger swings between them. 

We now turn to gauging the ability of our model to replicate various features of fi- 
nancial markets. One crucial aspect of financial markets returns is their volatility, as 
well as how this volatility evolves over time. Volatility is more relevant today than ever, 
with large spikes possibly occurring in short periods of time. Financial markets returns 
generally display volatility 'clusters'. These clusters indicate that once the volatility is 
high, it tends to remain high for a while, and that similarly, once it has come down, 
it tends to remain low for some time. A convenient and well-accepted way of modeling 
such characteristics is through the use of the Auto Regressive Conditional Heteroskedas- 
ticity (ARCH) model pioneered by Engle Q or through the use of the more encompass- 
ing Generalized Regressive Conditional Heteroskedasticity (GARCH) model proposed by 
BoUerslev [1]. 

Whether working on pricing a derivative product, attempting to hedge an exposure, 
optimizing a portfolio in a mean-variance framework, or estimating the Value- At-Risk of 
a position, the ability to capture and model the stochasticity and the clustering prop- 
erties of the volatility is paramount. Not doing so can lead to the wrong probability 
distribution being used, since volatility clusters impact the shape of returns distributions 
in two important ways. First, the fact that a period of calm statistically tends to be 
followed by another period of calm indicates that there will be a fairly large amount of 
probability mass around the mean (return) . Graphically this phenomenon translates into 
a probability distribution function that is higher than the Gaussian one in the vicinity 
of the mean. Second, the fact that a period of extreme movements statistically tends to 
be followed by another period of extreme movements indicates that there will be non- 
negligible amounts of probability mass in the tails areas. Graphically this phenomenon 
translates into a probability distribution function that is higher than the Gaussian one 
in the vicinity of the tails. 

The fairly wide-ranging GARCH specification models the volatility as both a function 
of past squared return shocks and of past levels of itself. If both past return shocks and 
past volatility levels are low, for instance, the odds are that the next volatility levels will 
remain low. If past volatility levels are low but recent past return shocks are high, the 
levels of volatility will likely increase. If past volatility levels are high but recent past 
return shocks arc low, the levels of volatility will perhaps decrease. Finally, if both past 
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return shocks and past volatility levels are high, the odds are that the next volatility 
levels will remain high. In a GARCH(p,q) specification, the number of lags p and q 
allowed for past shocks and past volatility levels are limitless. However, Hansen and 
Lunde [TtI] explore the ability of 330 different ARCH/GARCH models to capture the 
features of various financial returns and come to the conclusion that a GARCH(1,1) 
model performs just as well as the more 'sophisticated' ones. Our specification for the 
volatility of at time t is thus a GARCH(1,1) described by the following equation: 

(T^ =ao + aie?_i+/?i(T2_i, (28) 

where and <j1__i are the one-period lagged squared return shock and one-period 
lagged variance, respectively. 

We simulate 100 lattice time series for two different levels of our tuning parameter 
X- In the first set of simulations, x is set equal to a value of 0.0005 while in our second 
set of simulations, x is set equal to a value of 0.0013. This allows for the generation of 
different returns distributions and times series, namely financial markets returns with 
varying degrees of activity: The set of simulations using a x value of 0.0005 reflect a 
somewhat more volatile market than the set of simulations using a x value of 0.0013. 
Previous studies have shown that even a simple lattice model (9| is able to produce returns 
volatility dynamics displaying some ARCH/GARCH effects, and also that it is able to 
produce returns volatility dynamics that are rather consistent with those of NASDAQ 
historical returns. With the present model, we are here conducting a more extensive 
investigation: We estimate GARCH(1,1) fits from a large number of randomly chosen 
lattice configurations for the purpose of demonstrating the stability, consistency, and 
flexibility of our lattice model. 

The resulting parameters obtained are shown in Fig. [TUl and in Fig.[Tl]for x — 0.0005 
and X = 0.0013 respectively, for a total of 200 different lattice time series. For visual ease, 
the parameters are ordered before being plotted. In Fig. [TO] one can thus immediately 
see that the fii parameter is between 0.85 and 1.0 in about 95% of the cases, a result 
consistent with Nasdaq GARCH(l,l)-estimated figures in |9|. In Fig.[TT]the /?i parameter 
is between 0.85 and 1.0 in about 80% of the cases. This indicates that the lagged volatility 
feedback is strongly consistent across various market conditions and various simulations. 
In Fig. [TO] one can also see that the ai parameter is below 0.15 in about 99% of the 
cases, a result also consistent with Nasdaq GARCH(l,l)-estimated figures in [9[. Finally, 
in Fig. [TT] the ai parameter is below 0.15 in about 99% of the cases, although a higher 
proportion is above 0.05 when compared to the values it takes in Fig. [TO] These results 
indicate that the lagged return shock feedback is also strongly consistent across various 
market conditions and various simulations. The ao parameter is consistently very low in 
both cases, as it should at these high frequencies. 

It is also interesting to note that for each parameter, we obtain some 'outliers' in 
about 5% to 20% of the cases, although they are outliers only in the sense that they differ 
from the other 95% to 80% of the estimates that are themselves incredibly consistent, 
and are not outliers in the sense that their values would be considered too extreme or 
unreasonable. It is important to note that when running multiple simulations, one is 
bound to obtain some results that differ somewhat from the estimates' consensus. For 
instance, even if one was to simulate an exact GARCH(1,1) process many times and 
subsequently estimate the parameters from the simulated data, some percentage of the 
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Figure 10: GARCH(1,1) fit parameters for a sample of 100 lattice time series at x = 0.0005. The 
configuration counter c has been subject to sorting with respect to /3i. 
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Figure 11: GARCH(1,1) fit parameters for a sample of 100 lattice time series at x = 0.0013. The 
configuration counter c has been subject to sorting with respect to /3i. 
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estimated parameters would stray from the median estimates. In summary, 100% of 
our estimated parameters are sensible, across various market conditions (more or less 
volatile) and thus across varying returns distributions, and over a myriad of simulations. 

5. Conclusion 

The subject of this work has been to explore a class of models designed to simulate 
the properties of financial markets. The output of the model is a time series of returns, 
from which gains distributions and related features could be derived. The potential of 
the model for replicating market dynamics, as described by standard financial analysis 
tools, was the primary aim of this study. 

The production of our time series has been done by numerical simulation based on a 
lattice description of fields in time-asset space. We have restricted ourselves to a one-asset 
model linked to an interest rate account. The dynamics are based on a gauge invariant 
lattice action which, when quantized, gives rise to eliminating arbitrage opportunities up 
to stochastic fluctuations, thus reflecting real market conditions. The second pillar of 
the model is an updating prescription that evolves the lattice fields into a self-organizing 
critical state. This appears to be an essential element for reproducing certain stylized 
features of real markets. 

As a third feature, a parameter has been introduced as a tuning tool, through which 
a variety of market characteristics, ranging from quiescent to volatile markets, can be 
modeled. 

An extensive analysis of a very large number of time series features evaluated by 
a GARCH(1,1) analysis was performed. It turns out that close to 100% of the lattice 
model-generated time series give rise to sensible analysis parameters, rendering the model 
results almost indistinguishable from historical market data. In particular we could verify 
this observation across various market conditions and varying returns distributions. 

We conclude that the model shows promise as a modeling tool for financial time series 
and look forward to further development and applications. 

Appendix A. Gauge fixing 

In the notation of 8J, the probability density for a particular gauge field component 
has the form 



where Lq and Lq are positive coefficients independent of 6*^ (x) . They reflect the (local) 
environment of the link variable. Under a gauge transformation, writing 



Pe{&n{x)) oc exp(-/3(Le exp(6'^(x)) + cxp(-6'^(a;))ie) ) 



(A.l) 



g{x) = e'^(-) , 



(A.2) 



we have 



(A.3) 
(A.4) 
(A.5) 
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The transformation laws for Lq and Lq can be derived directly by an examination of 
the lattice action S'[0, $]. They are also obvious from the fact that the action is gauge 
invariant and the arguments of the exponential functions of (|A.1[) are made up from 
invariant contributions to it. We now choose the gauge transformation by requiring the 
new coefficients to be equal 

^ ^ ^2h(.)Le^-2h(.+e,) ^ ^ ^ (A.6) 



or 



A solution of (IA.7I) is 



hix + e,)-h{x) = llogi^). (A.7) 



h{x + e,) = A±liog(||) (A.8) 

h{x) = ^log(||), (A.9) 

where A is a real parametefl On all other sites, besides x and x + e^, the gauge 
transformation function hix) is arbitrary. (For definiteness one may choose it to be 
zero.) In the new gauge, using (jA.6|) . we have 

+ e-<^»L'Q ^ 2yi^cosh(f?;(x)) . (A.IO) 
Thus, the probability distribution function is 



p's{e[Xx)) « exp(-2/3^L'eZ^cosh(0;(a;)) ) . (A.ll) 

Dropping the primes gives ((T5| . Citing LqLq = LqL'q, we note that the variance of the 
the probability distribution function is not altered by the gauge transformation. 

Again, in the notation of [SJ, the probability density for a particular matter field 
component has the form 

p^{(j>^{x)) oc exp(-/3(L$ exp((/)^(a::)) + cxp(-(/)p(a;))i$) ) (A. 12) 

where L$ and Z$ are positive coefHcients independent of 0^(x), reflecting the (local) 
environment of the field variable. In this case, changing the gauge (IA.2I) entails the 
transformations 

(/)(a;) ^ ^'{x) = h{x) +4>{x) (A.13) 
L'^ = e''(^)L$ (A.14) 
Li L'i = Z<[.e-''(^) . (A.15) 



^It may be used to control the effect of the gauge transformation on the matter fields (p{x) and 
4>(x + e^). 
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The requirement L'^/L'^ = 1 leads to 

h{x)^-\\og{^), (A.16) 

while the gauge function is arbitrary on all sites other than x. Proceeding in the manner 
above we have 

L'^e'^''^^^^ + e-^^(-)L'^ = 2^IjJ7^cosW^{x)) , (A.17) 
and thus obtain the probability distribution function for the matter field 

Mc^'^ix)) cx eM-'^P^ L'^L'^cos\v{c^'^{x)) ) . (A.18) 

Dropping the primes gives p^ . Again, because of L^pLip — L'^L'^, the gauge transfor- 
mation does not change the variance of the distribution. 
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